Theory of Cherenkov radiation in periodic dielectric media: Emission spectrum 
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The Cherenkov radiation is substantially modified in the presence of a medium with a nontrivial 
dispersion relation. We consider Cherenkov emission spectra of a point or line charge, respectively, 
moving in general, three- (3D) and two-dimensional (2D) photonic crystals. Exact analytical ex- 
pressions for the spectral distribution of the radiated power are obtained in terms of the Bloch 
mode expansion. The resulting expression reduces to a simple contour integral (3D case) or a one- 
dimensional sum (2D case) over a small fraction of the reciprocal space, which is defined by the 
generalized Cherenkov condition. We apply our method to a specific case of a line source moving 
with different velocities in a 2D square-lattice photonic crystal. Our method demonstrates a rea- 
sonable agreement with numerically rigorous finite-difference time-domain calculations while being 
less demanding on computational resources. 
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I. INTRODUCTION 

Back in 1934 Cherenkov reported the observation of 
the electromagnetic radiation produced by an electron 
moving in a dielectric medium at a velocity greater than 
the phase velocity of light in this medium Such a 
radiation possesses a unique angular and frequency spec- 
trum and is called Cherenkov radiation [2j]. A nontrivial 
dispersion relation of a medium leads to substantial mod- 
ifications of the Cherenkov radiation. It has been shown 
that an electron moving in a homogeneous medium with 
dispersion should emit at any velocity 0|. Richer spatial 
distribution of the emitted radiation including intensity 
oscillations behind the Cherenkov cone is a signature of 
the radiation in such a medium [1, [H, [f| . 

To understand the properties of the Cherenkov radia- 
tion one can represent the moving electron with space- 
time dependence of the corresponding current density 
J (r, t) ~ 5 3 (r — vt) as a superposition of plane waves 
S 3 (r — vt) = J2k ex P (*k ■ r — ik ■ vt) with different wave 
vectors k and frequency k • v, where v is the electron ve- 
locity. Only plane waves with frequency and wave vector 
fitting the medium dispersion uj (k) can resonantly ex- 
cite electromagnetic modes in the medium, which gives 
the Cherenkov resonance condition Q: 



u> (k) = k • v. 



(1) 



In a homogeneous, non-dispersive medium with refrac- 
tive index n, the dispersion relation is simply given by 
uj (k) = (c/n) |k| and the Cherenkov condition (J) leads 
to a well known conical wave front with an aperture 
cos0 = c/ (n |yj) and a condition on the electron veloc- 
ity |v| > c/n [2 0|) c being the vacuum speed of light. 
In an inhomogeneous medium the interplay between in- 
terference and propagation can result in an engineered 
nontrivial dispersion relation uj (k). For example, peri- 
odic dielectric media (photonic crystals) @, S| substan- 
tially modify both dispersion and diffraction of electro- 
magnetic waves possessing many unusual and novel op- 



tical phenomena, including modification of emission dy- 
namics pjl El El, ultra- refraction [H, [3, [H, EDI and 
photon focusing [l7l. ll8L effects. The present work 
focuses on an analytical understanding of the influence 
of a periodic medium on the Cherenkov effect. 

Several studies on the modification of the radiation 
produced by a charged particle moving near or inside 
periodic dielectric media are available. The Cherenkov 
radiation in cholesteric liquid crystals has been analysed 
in Ref. [2(J. The modification of the Smith-Purcell ra- 
diation has been recently studied both theoretically and 
experimentally near a surface of a two- (2D ) a nd three- 
dimensional (3D) photonic crystal [gj Hi (H H3]. 
The Cherenkov radiation generated by an electron mov- 
ing inside an air pore of a 2D photonic crystal perpen- 
dicular to the periodicity plane has been used to map 
its photonic band structure in Refs. [26|, [27]]- In all 
above mentioned reports, the theoretical analysis of the 
Cherenkov effect has been done in the plane wave ba- 
sis. Spatial and spectral modifications of the Cherenkov 
radiation produced by an electron moving in the period- 
icity plane of a 2D photonic crystal have been studied in 
Ref. [23| using the finite-difference time-domain (FDTD) 
method. To date, there do not exist any reports on the 
general theory of the Cherenkov effect in an arbitrary 3D 
periodic dielectric medium [3^ |. 
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Figure 1: A sketch of a periodic medium and a charge trajec- 
tory. Basis vectors of the lattice are shown. The coordinate 
system is chosen with one axis along (fy) and the other per- 
pendicular (fx) to the charge trajectory. 
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The main purpose of the present work is to develop 
such a theory and to provide a simple expression for the 
Cherenkov emission spectrum (energy loss spectrum) for 
a point or line charge, respectively, moving with velocity 
v in an arbitrary direction inside a general, 3D or 2D 
photonic crystal. To achieve this goal we derive an ana- 
lytical expression for the power emitted per unit length 
of the charge trajectory in terms of the Bloch mode ex- 
pansion. This expression is further reduced to a simple 
contour integral (3D case) or a one-dimensional sum (2D 
case) over a small fraction of the reciprocal space. As 
a result, to calculate the Cherenkov emission spectrum, 
Bloch eigenmodes and their corresponding group veloci- 
ties are required only along an integration path (3D case) 
or at a discrete set of Appoints (2D case), considerably 
reducing computational demands. The integration path 
and the discrete set of points are defined by the general- 
ized Cherenkov condition. Our theory confirms that the 
Cherenkov radiation does exist in a periodic medium for 
an arbitrary electron velocity It also predicts an 

enhancement of the radiated power near the frequencies 
corresponding to the vanishing component of the group 
velocity, which is orthogonal to the electron trajectory. 

The paper is organized as follows. In Section [II] the 
general solution of Maxwell's equations is summarized 
for an arbitrary periodic medium. In Section IIIII an an- 
alytical expression is derived for the power radiated per 
unit length by a moving point charge in 3D and a line 
charge in 2D periodic media. In Section [TV] we apply our 
theory to calculate the Cherenkov emission spectra in the 
particular case of a 2D photonic crystal. The predictions 
of the analytical theory are substantiated by numerically 



rigorous FDTD calculations. Section [V] concludes the 
paper. 



II. RADIATED FIELD 

We consider a point (line) charge q uniformly moving 
with a velocity v along some direction in a general, infi- 
nite periodic 3D or 2D medium e (r) = e (r + R) (Fig. 
[J). Here R is a vector of the direct Bravais lattice, 
R = J2ili a ii h ls an integer and en is a basis vector 
of the periodic lattice. It is assumed that the medium is 
linear, nonmagnetic and that no absorption takes place. 
Then the relevant Maxwell's equations read in SI units: 

VxE(M) = - M o^H(r,t), (2) 
V x H (r, t) = e Q e (r) J^E (r, t) + J (r, t) , (3) 

where, the electric (magnetic) field is denoted by E (H). 
An electromagnetic field is produced by a current source 
J(r,i), which in the case of the moving point (line) 
charge is defined as 

J (r, t) = qvS d (r - vi) . (4) 
Here d — 2, 3 is the dimensionality of the periodic lattice. 
In the frequency domain, a general solution of Maxwell's 
equations ((2][3j) for an arbitrary current source J (r, t) and 
a periodic dielectric function e (r) is given in terms of the 
Bloch eigenmode expansion 0, l3(i| 
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Here J (r, oj) is the Fourier transform of the current den- 
sity J (r,t), for the moving point (line) charge Q given 
byfia 

J(r,w) = qf||(5 !i (r_ L )exp f iw rf|) • ( 6 ) 

The coordinate system is chosen with one axis, fy, be- 
ing parallel and other, orthogonal to the elec- 

/^N (LA 

tron trajectory (Fig. [J). E^ n (r) and E^J (r) are gen- 
eralized transverse and longitudinal Bloch eigenmodes 



fgl. f30l| characterized by the band index n, the wave vector 
k and the eigenfrequencies uj^ and w^n , respectively. 
The Bloch eigenmodes satisfy standard lattice periodic 
boundary conditions. As it will be shown in the next 
Section, only transverse Bloch eigenmodes contribute to 
the Cherenkov radiation field. The asterisk (*) and <g> de- 
note the complex conjugate and the outer tensor product 
in 3D space, respectively. Bloch eigenmodes satisfy the 
homogeneous wave equation and fulfill the normalization 
conditions 



J d d re(r)E^(r)-E 
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and completeness relations 



[ d d k vMrMrOE^ (r) ® E& } * (r') = (27r) d i<5 d (r - r') , 



(«) 



r 



where a, /3 = T or L, and 1 is the unit tensor. In Eqs. ([5]) with a total power radiated per frequency interval 
the /c-space integration is performed over the first Bril- [oj, uj + dui] given by 
louin zone (BZ) of the periodic medium and the sum- 
mation is carried out over different photonic bands. A 
positive infinitesimal 7 in J5]) assures causality [§]. 



III. EMISSION SPECTRUM 



-Re 

7T 



Vo 



(12) 



The emitted power of the Cherenkov radiation in a 
general periodic medium is given by the rate at which 
the moving charge does work on the surrounding electro- 
magnetic field. For an arbitrary current density J (r, t) 
in a 3D or 2D volume Vq, the time-dependent emitted 
power is given by [31| 

P(t) = - J d d rJ(r,t)-E(r,t). (9) 

Vo 

The total energy U radiated by the current J (r, t) is 
obtained by integrating |(9]) over all moments of time 

/oo 
dtP(t). (10) 
-00 

The time integral in (fTU|) can be further transformed into 
the integral over frequency (see Appendix A) 

U= dujP(uj), (11) 
Jo 

I 



To obtain the power emitted per unit length of the charge 
trajectory the integration volume Vq should be chosen as 
a cylinder coaxial with the trajectory, while the integral 
itself should be normalized by the cylinder length I. In 
the 2D case, the volume integral is reduced to the surface 
integral over a rectangle coaxial with the charge trajec- 
tory and the result should be normalized to the rectangle 
length. 

We further derive the spectral dependence of the power 
(dP/dl) lfT2]) radiated per unit length by the point (line) 
charge |(4]) uniformly moving in a periodic medium. As- 
suming that the presence of the moving charge does 
not change the band structure of the periodic medium, 
the electromagnetic field E (r, cj) surrounding the moving 
charge can be expressed in the form of the Bloch eigen- 
mode expansion ((5j) . This expansion is valid for any point 
r in the medium being different from, but as close as re- 
quired to, the charge trajectory. Substituting the Fourier 
transform of the current density |(6|) and the Bloch mode 
expansion (0 in the equation (fT2"| we obtain 



dP 
~dJ 
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(13) 
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where a = T, L. We have readily performed the 
space integration in the transverse direction f j_ and used 
the Bloch theorem E^ (r) = (r) exp (ik • r) , where 

e kn ( r ) i s a lattice periodic function. To avoid having to 
deal with the "bremsstrahlung" radiation we limit our- 



selves to the electron trajectories which do not cut dielec- 
tric interfaces in the periodic medium. Such trajectories 
are necessarily rationally oriented with respect to the pe- 
riodic lattice. In this case the function (e^ (r\\) ■ in 
(fl4| as well as its complex conjugate are both one di- 
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mensional periodic functions with a period a defined by 
a particular orientation of the electron trajectory. Then, 
Eq. (fl4|) can be further simplified to (see Appendix B) 



= 27rg 2 



,(«) 
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(15) 

Here fc|| is the component of the wave vector paral- 
lel to the electron trajectory. c TO (k;n) is the m-th 
(to S Z) Fourier coefficient of the periodic function 



e kn ( r " I ' *n I defined as 



^ n ) = - a [ dr W { e ^(r\\)-r\i)e-^K (16) 



Taking into account the expression (fl5"l) and the rela- 
tion Re [i (Re [z] + i Im [z])] = — Im [z], the power radiated 
by a moving charge per unit length is given by 



dZ (2 7r ) d - 1 7re isz V " |v| a 



Im 



(k;n) 

This expression can be further integrated along the direction k\\ in the A:-space yielding 
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(17) 



(18) 



In the 3D resulting surface integral is taken over 

a plane S. In the 2D case, the integration reduces to 
an integral over a line C (Fig. [2]) . Both the integration 
plane and the integration line should be orthogonal to 
the electron trajectory and are defined by the following 
relation 



fc|| = — + —TO. 

" v a 



(19) 



Here the integer to should be chosen in such a way that 



the wave vector ku stays in the first BZ. Further, taking 
the limit 7^0+ and using the relation 



Im 



lim 7=7- 



17 



-7rf(w±w£?), (20) 



the spectral radiated power l|18p can be expressed in the 
form 



dP 



Luq 



(27T) 
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E 



2cj 



(T) 
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,(T) 



(k; n)| Z (5 ( w - 4^) - 5 (a, + + ^ |4 L) (k; n) [ «S ( W ) 



(21) 

The eigenfrequencies of the Bloch modes are positive [8J|, so the second term in Eq. lj2l"j) containing the delta 



function 6 



is zero for all frequencies. The third term in Eq. (J2TJ) is due to the work the current does on 

the longitudinal part of the electromagnetic field. In the presence of free charges the longitudinal part of the field 
corresponds to the static electric field and the work done against it results in non-radiative energy transfer with a 
nonzero contribution only at zero frequency. In what follows we will disregard this non-radiative contribution and 
will limit ourselves to the radiation into propagating electromagnetic waves only. Then the spectral radiated power 
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is given by 



dP_ 
~df 
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Figure 2: (Top) Diagram to define the integration plane 5 and 
the integration contour OS (dashed line) in Eq. (| 1 8 124 j) . Iso- 
frequency surface enclosed in the first BZ of the FCC lattice 
is shown for a normalized frequency Ukn = u> inside the first 
bandgap of a 3D inverted opal [321 ] . (Bottom) Diagram to 
define the integration line C and the set of points {&x,i} (two 
thick dots) in Eqs. I|18l28p . Iso-frequency contour enclosed in 
the first BZ of a square lattice PhC is shown for a normalized 
frequency ojkn = uj inside the first bandgap. The plane S and 
the line C are defined by the relation fey = 5k\\ = j^j + fi. 
The choice of the coordinate system with one axis, fty , parallel 
to the electron trajectory is shown. 



The argument of the Dirac delta function in Eq. lf22|) 
is a function of the wave vector. One can use this fact to 
further reduce the dimensionality of the (d — 1) fc-space 
integral. In the 3D case, using the relation 



/ 



d d kf(k)S(g(k)) 



jd-li 



av 
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where dV is (d— 1) dimensional surface defined by 
,g (k) = 0, the integral over the plane S is converted into 
a contour integral 
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The contour dS is defined by the relation (fl9f and 



(24) 



(25) 



It is an intersection of the iso-frequency surface with 
plane S (Fig. El-top). In the 2D case, the relation 



*(/(*)) = £ 



<S(fc- A;«) 
!/'(**)! 



(26) 



can be used, where summation is taken over all solutions 
of the equation / (k) = 0. Substituting (|26|) into ([22)1 we 
obtain 



2D 



47T£( 



Cm' (k;ra) 




- kj_,i) 




aw, 


2/n± 


) 


fc_L,i 



(27) 

where {k±,i} are simultaneous solutions of the equa- 
tions (fl9|) and f25|) given by the intersections of the iso- 
frequency contour with the line C (Fig. [2]-bottom). Per- 
forming fc-space integration, we finally obtain 



dP 



21) 



47T£( 



E 



S T ) 



(k;n) 



(28) 



where the function in brackets is calculated for the wave 
vectors corresponding to the set {fcj_,i}- 

Formulas l(24|) and ([28)1 constitute the main result of 
the present work. They give the power radiated by the 
moving point charge (3D) or line charge (2D) q in the 
spectral interval [u, u + duj] per unit length of the tra- 
jectory for a 3D and 2D periodic medium, respectively. 
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The radiated power is proportional to the Fourier co- 
efficients e„ (k;n), which effectively describe the local 
coupling strength between the current density produced 
by a moving charge and the electromagnetic field at the 
electron location. The gradient and derivative of the dis- 
persion relation v 9 = Vk ± w^ and v 9 ^ = duj^/dk± 
yield the component of the group velocity, v 9 , of the 
Bloch eigenmode (k;n), which is orthogonal to the elec- 
tron trajectory. The Cherenkov radiated power is pro- 
portional to the inverse of this component of the group 
velocity. That means that the radiated power can be 
strongly enhanced not only if the group velocity itself 
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Figure 3: Diagram to illustrate the generalized Cherenkov 
condition l|29ll . A 3D representation of the photonic band 
structure (top) of the 2D PhC (bottom) is shown for TE 
polarization. An infinite 2D square lattice of air holes in a 
dielectric medium is considered. The radius of the holes is 
r = 0.4a, the dielectric constant of the background medium is 
£ = 12.0. Only the first band in the first BZ is presented. The 
right-hand side of Eq. lf29ll defines the set of planes for dif- 
ferent m. The intersection (dashed line) of these planes with 
the band structure determines the Bloch modes contribut- 
ing to the Cherenkov radiation. Here it is supposed that a 
line charge moves along the 2-axis in the crystal with velocity 
Ivl = 0.15c. 



is small for some frequency, but also if the component 
of the group velocity orthogonal to the electron trajec- 
tory becomes small. At the same time suppression of the 
Cherenkov radiation is possible if for some frequency the 
current density produced by a moving charge is not cou- 
pled to the corresponding Bloch mode and the Fourier 
coefficients c m (k; n) is small. 

Only eigenmodes with the wave vectors on the contour 
dS (|24|) and from the set |28|) contribute to the 

radiated power at a given frequency. It is important to 
realize that Eqs. lfl9|) and (|25)l defining the contour dS 
and the set {fcj_,i} are equivalent to the Cherenkov reso- 
nance condition (TJ. In fact, substituting ([25)) in (|T9j) and 
taking into account that the scalar product in {TJ results 
in v • k = |v| A; 1 1 one obtains the generalized Cherenkov 
condition for a periodic medium 



CO 



(T) 
kn 



Ivl fell 



|v| — m. 



(29) 



In the 4D (3D) (w-fe)-space the right-hand side of the 
relation (|29l) defines a hyperplane (plane) whose inter- 

(T) 

section with the band structure, , determines Bloch 
modes contributing to the Cherenkov radiation (Fig. [31 
top). Nonzero integers m ensure that such an intersec- 
tion and consequently the Cherenkov radiation exist in a 
periodic medium for an arbitrarily small charge velocity. 
In a homogeneous medium m — and the Cherenkov 
condition reduces to a standard form Wk = |v| fen. 

As a simple check of our theory we show in the follow- 
ing that the final formulas <|24I28|) reproduce the limit of a 
homogeneous medium with the dielectric constant e. For 
a given frequency to, the wave vector |k| and the group 
velocity |vj_ | are given by |k| = (w^/e) /c and |vj_| = 
(c|kj_|) / (y/s\\s.\) respectively, with k± = k k|| being 
the component of the wave vector perpendicular to the 
electron trajectory. The appropriately normalized eigen- 
modes are plane waves E = (l/^fe) eexp (ik ■ r), where e 
is a polarization unit vector orthogonal to the wave vec- 
tor k. Further, according to the Eq. (fl9|) the wave vector 
component fey is equal to kn = to/ |v| with m — and the 
coefficient cq is given by cq = |kj_| / (y/e |k|). Then in the 
3D case, taking into account that an integration contour 
dS is a circle of radius |kj_| and performing integration in 
polar coordinates with dk = \k± \ d<j>, the radiated power 
(J24J) is given by 
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(30) 



which finally yields the usual results of the Frank- Tamm 
theory Q 



dP 
~dJ 
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(31) 
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In the 2D case Eq. J28]) yields 

dP\ 2D _ _L <f_ 

dl J h 2ne c^/s 1 



e v 



(32) 



IV. NUMERICAL RESULTS 

In this section the analytical results developed in the 
previous section are applied to the numerical study of the 
Cherenkov radiation in a 2D photonic crystal. An infinite 
2D square lattice of air holes in a dielectric medium is 
considered. The radius of the holes is r = 0.4a, while the 
dielectric constant of the background medium is e — 12.0. 
A line charge oriented perpendicular to the periodicity 
plane of the crystal moves along the z-axis with a velocity 
v, staying always in the space between air holes (Fig. [3] 
bottom). The corresponding current density, Eqs. (@]) 
and ((6]), generates an electric field ((5j) polarized in the 
periodicity plane (transverse electric or TE polarization) 
[1,0], i.e., the Bloch eigenmode expansion should include 
the TE polarized Bloch modes only. The first TE band 
for the considered PhC is presented in the figure [3]-top. 
The band structure was calculated using the plane wave 
expansion method [33]. 

To find the power radiated by a charge moving with a 
given velocity v, all Bloch modes contributing to the ra- 
diation should be determined. These modes are specified 
by the solutions of the relation (j2"9]l . In what follows we 
restrict our analysis to the frequency range of the first 
band of the considered PhC structure. In figure [3]-top, 
solutions of the Cherenkov relations (|29|) are graphically 
illustrated for m = 0, —1, —2 (dashed lines) and charge 
velocity v = 0.15c. The frequencies satisfying relations 
(|29|) determine the spectral range of nonzero contribu- 
tion to the Cerenkov radiated power, the Cherenkov band. 
The evolution of the Cherenkov band is presented in fig- 
ure S] as a function of the charge velocity. 

For charge velocity v = 0.15c, the Cherenkov spectrum 
is given by the intersections of the band structure with 
the planes m = — 1 and m = — 2. The plane correspond- 
ing to m = does not intersect the band structure of the 
crystal (Figs. [3j [4]). For smaller charge velocities, more 
and more planes intersect the photonic band structure, 
and the Cherenkov spectrum is built from a number of 
discrete sub-bands. For sufficiently small charge veloci- 
ties the spectral range of the first photonic band becomes 
densely filled with the discrete sub-bands (Fig. [4j) . 

In the long wavelength limit the periodic medium is 
effectively homogeneous. For the PhC considered in 
Fig. H] the effective refractive index is equal to n c s = 
^/ScB ~ 2.186. Consequently, for m — the relation 
(|29ll imposes a condition on the minimal charge veloc- 
ity to produce Cherenkov radiation at small frequencies, 
namely v > v m - m — c/n c g 0.457c. For charge veloci- 
ties larger than the threshold value the Cherenkov band 
covers the spectral range from zero to the maximum fre- 
quency, which is defined by the intersection of the band 
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Figure 4: Cherenkov radiation band. Sub-bands denned by 
the intersections of the band structure with the planes cor- 
responding to the different m's are shaded in light and dark 
gray. Vertical lines mark the charge velocities used in the fur- 
ther calculations, v = 0.15c, v — 0.3c and v = 0.6c in the top 
panel and v = 0.1c in the bottom panels, respectively. 



structure with the plane m = at the first BZ boundary 

(Fig.rj. 

To compute the radiated power from Eq. (f28| . one 
should calculate the Bloch modes along a charge trajec- 
tory, their Fourier transforms and corresponding group 
velocities for wave vectors belonging to the intersections 
defined by Eq. I|29p . The calculation of the Cherenkov 
spectrum is illustrated in figure [H for a line charge (q = 
1.6 x 10~ 19 C) moving with the velocity v = 0.1c. To cal- 
culate Bloch modes and group velocities, the plane wave 
expansion method [33| and the Hellmann-Feynman theo- 
rem were used, respectively. For the velocity v = 0.1c the 
Cherenkov spectrum consists of three sub-bands defined 
by the planes m = —1, m= —2 and m = —3 (FigsH] [oj, 
respectively. The Fourier coefficients, Cm (k;n), (Fig.[5]- 
b) and the orthogonal component of the group velocity, 
v 9 ± , (Fig. [5]-c) are nonzero only within the sub-bands. 
Both Fourier coefficients and the orthogonal component 
of the group velocity approach zero at sub-band edges 
A, B and C. At the edges D, E, and F only the orthog- 
onal component of the group velocity is zero, while the 
Fourier coefficients have finite nonzero value. Calculation 
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Figure 5: Cherenkov emission spectrum for the charge veloc- 
ity v — 0.1c. Projections of the first band of the considered 
photonic crystal and the planes m = — 1, m = —2 and m = —3 
defining the Cherenkov band on the k x -ui plane are shown 
(a). The Fourier coefficients, Cm (k;n), the orthogonal com- 
ponent of the group velocity, v 9 ± , and the Cherenkov power 
spectrum are shown in panels (b), (c) and (d), respectively. 
Vertical lines mark the edges of the Cherenkov sub-bands. 
Contribution from the sub-bands corresponding to m = — 1, 
m = —2 and m = —3 are shown as dashed, solid and dashed- 
dotted lines, respectively. 



of the Cherenkov radiated power at the band edges A, B 
and C leads to the indeterminate limits of the form 0/0, 
which can be evaluated using l'Hopital's rule and is equal 
to zero. At the band edges D, E and F the Cherenkov 
power diverges in an intgegrable way. 

In figured] the Cherenkov radiated power spectra are 
shown for charge velocities v = 0.15c, v = 0.3c and 
v = 0.6c. For charge velocities smaller than the threshold 
value v m i n w 0.457c the Cherenkov radiation is nonzero 
only within single or multiple spectral bands. For veloc- 
ities above the threshold, the radiated power is nonzero 
almost everywhere within the first band, approaching 
asymptotically the value of the Cherenkov radiated power 
in a homogeneous medium with n — n e g for small fre- 
quencies. The radiated power calculated using Eq. f32|) 
for v = 0.6c and n e s — 2.186 is shown in figure [6] (bottom 
panel) as dotted line. The Cherenkov radiated power is 
enhanced near those frequencies where the group velocity 
component orthogonal to the charge trajectory vanishes, 
while the Fourier coefficients remain finite (Fig. [6]) . 

To substantiate our analytical results the direct nu- 
merical integration of the Maxwell's equations has been 
performed using rigorous finite-difference time-domain 
(FDTD) method [3J|. The simulated structure was a 
10a x Na lattice of air holes in a homogeneous medium 
with e = 12.0. The longitudinal size of the periodic struc- 
ture was set to N = 188, N = 376 and N = 752 lat- 
tice constants for an charge velocity v — 0.15c, v = 0.3c 
and v = 0.6c, respectively. The lattice was surrounded 
by a 2a wide layer of homogeneous material. The sim- 
ulation domain was discretized into squares with a side 
A = a/18 and was surrounded by a 35-cell-wide perfectly 
matched layer (PML) [35| . The time step of integration 
was set to 98% of the Courant value. The moving line 
source (@]) was modeled as a current density source j34[ 
with the Dirac delta function represented via an appro- 
priately normalized Kronecker delta Sij/A . The charge 
trajectory was oriented in the longitudinal direction of 
the system, placed in the geometrical center of the crys- 
tal, exactly between the 5th and the 6th row of holes. 

For the numerical computation of the radiated power 
the electric and magnetic fields were stored at a detec- 
tor surface enclosing the crystal, and their Fourier trans- 
forms with respect to time were found by discrete Fourier 
transformation. The longitudinal dimension of the struc- 
ture was different for different charge velocities in order 
to keep the integration time at the detector and conse- 
quently the spectral resolution constant. The detector 
surface was situated in the close vicinity of the crystal 
boundary. The total radiated power per unit length was 
then calculated as 



dP 



1 2 
d 



LJ 

— I dz S (z, to) ■ n 

71" J 



(33) 



where S (z, to) = ^Re [E (z, oS) x H* (z, to)] is the Poynt- 
ing vector at radiation frequency u, D is the length of 
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Figure 6: Cherenkov emission spectra for different charge ve- 
locities in a 2D photonic crystal. Solid lines correspond to 
the analytical results (l28l) . Dotted line in the bottom panel 
(v = 0.6c) corresponds to the Cherenkov radiated power 
in the homogeneous medium with n c s = 2.186, Eq. lf32]l . 
Radiated power spectra obtained using FDTD method are 
shown as dashed lines. For charge velocity v — 0.15c radiated 
power spectra are show for both 10a x 188a (dashed line) and 
20a x 188a (dashed-dotted line) structures. 



the detector surface along the trajectory and n is a unit 
vector orthogonal to the detector interface. 

An overall very good agreement between the results of 
the analytical (Fig. [6j solid lines) and numerical calcula- 
tions (Fig.El dashed lines) is obtained. Both the spectral 
range of a nonzero radiated power and its absolute value 
are well represented using the FDTD method. The main 
difference can be traced back to Fabry-Perot-like oscilla- 
tions of the radiated power due to the finite-size effects 
in the FDTD calculations. In the finite structure, the 
Cherenkov radiated power stays considerably enhanced 
near the band edges having large but finite value. The 
total power oscillates around the analytical value becom- 
ing partially suppressed or enhanced for different frequen- 
cies. To confirm that these oscillations indeed result from 



the finite transverse dimension of the considered photonic 
crystal, we have performed simulations for the crystal 
with a double thickness (20a x 188a) for the charge ve- 
locity v = 0.15c. Resulting radiated power spectrum is 
shown in the corresponding panel in figure[6]as a dashed- 
dotted curve. One can see twice as many oscillations as 
in the case of the thinner structure, which is a typical 
signature of the Fabry-Perot like phenomena. The fur- 
ther enhancement of the radiated power in comparison to 
the infinite structure can be associated with the longer 
interaction time of the charge at resonance frequencies 
with the effectively slow Fabry-Perot modes of a pho- 
tonic crystal slab. 



CONCLUSIONS 



In this paper, analytical expression for the Cherenkov 
power emitted per unit length of the charge trajectory 
in the case of a general 3D and 2D periodic dielectric 
medium has been derived. The obtained formula for 
the Cherenkov power involves the calculations of Bloch 
modes and corresponding group velocities at a limited 
number of points of the reciprocal space only, making 
the application of the proposed method computationally 
not demanding. All calculations have been performed 
on a desktop PC and our method requires 5 to 10 times 
less CPU time then FDTD method. The analysis of the 
Cherenkov emission spectrum in the periodic medium 
reveals that the Cherenkov effect indeed exists for ev- 
ery electron velocity. Similar to the case of the modifi- 
cation of the dipole emission in a photonic crystal, the 
Cherenkov radiation can be suppressed if the coupling of 
the current density produced by a moving electron with 
a Bloch mode is poor. At the same time, an enhance- 
ment of the Cherenkov radiation is possible also if only 
the component of the group velocity orthogonal to the 
electron trajectory is small. We have illustrated the de- 
veloped analytical method and its conclusions using a nu- 
merically rigorous finite-difference time-domain method 
in a special case of a 2D photonic crystal and demon- 
strated a reasonable agreement between numerical and 
analytical results. 



Acknowledgments 



We are grateful to Sergei Zhukovsky for fruitful dis- 
cussions and careful reading of the manuscript. Finan- 
cial support from the Deutsche Forschungsgemeinschaft 
(DFG FOR 557) is gratefully acknowledged. 



10 



Appendix A 

Using the Fourier representation of the time depended 
real function F (r, t), 



F(r,i) = Re 



1 

2^ 



dco F (r, iS) e 



(34) 



= — / did (F (r, u) e- iut + F* (r, u) e luJt ) , 

J — oo 

for the electric field E (r, t) and the current density J (r, t) 
in Eq. ([9]), the total radiated energy (fTOf can be written 
in the form 



U 



r Re 



dt i duj 



(47T) 

{j(r,w) • E (r, 17) e- i(w+n) * + 
J (r, lj) ■ E* (r^e"*"-^*} 



dQ x 



(35) 



Changing the integration order and using the integral 
relation for the Dirac delta function 



1 r°° 

S(x) = — dye~^ 

27r J-oc 



we obtain for total radiated energy 
1 



(36) 



U 



47T 



-Re 



d' V / 

oo 



dfl x 



{J(r,w) ■ E (r, Q)5(w + Q) + 
J(r,w) -E* (r,0) *(&;-«)}] 



(37) 



Further, integrating over and using the symmetry of 
the Fourier transform of the electric field, E (r, — oj) = 
E* (r, w), the total radiated energy can be written as an 
integral over frequency (jTTJ) 



2 f°° 1 f 

U = / da; -Re / Aj(r,w)-E*(r,w) 

n Jo 2 J 

LVo 



(38) 



with coefficients c™' (k; n) defined in Eq. (fl6| , integrals 
lj and 4"' m (US can be transformed to 



(«) _ V n( 



oo 

«>(k;n) | 



dm e" 



and 



1/2 

oo ' 

4 Q) = E 4 Q) ( k ;™)7 / rf me i(fc|1 



i ( fc H~^~TT TO ) r ll 



(40) 



I p) r n 



(41) 



-1/2 



respectively. In Eq. (|40|) . integration over rn immediately 
yields 



J< a >=27r £ ^(kjnj^fcu-^r-^m). (42) 

TTl.— — (Yl \ 1 I / 



Integral in (HU is equal to I, if — t^t - = 0. Over- 
wise it results in 



< < fc 'l - H - TP (*/2) 



|v| 



(43) 



Then, in the limit I — > oo, Jj vanishes for — r^y — 

^-p ^ 0, while is equal to 2-k X^L-oc °p (^i n ) f° r fy| — 
t^t — ^-p = 0. Finally, using the function 



S(x) 



1, x = 

0, 2) ^ ' 



(44) 



The integrand in Eq. I|38|) coincides with the time- 
averaged radiated power of the monochromatic source 

J(r,w) [11. 



Appendix B 

Expanding a periodic function (r*||) • f||^ in the 

Fourier series 



oo 

(eS(r||)-fy)= E 4 Q) (k;n)e iS 



relation (fT4l) can be written in the form 



oo oo 



= 2tt g 2 E E cW*(k;n) C W(k;n)x (45) 

771= — GO J>= — OO 

<5 I fen — 7—7 — — m ) J (fcii - -p7 - — P 
v a / V 11 v a 



(39) which is non-zero only for m = p yielding relation (fTS)) . 
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